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ABSTRACT To investigate the effect of the dynamic deflection process of a two-dimensional vector 
nozzle on its actuator load and engine performance, a comprehensive modeling method of two-dimensional 
vector nozzle/actuator load/turbofan engine is proposed. Firstly, based on the two-dimensional vector 
nozzle spatial motion model, combined with force and moment analysis, the mechanical model of vector 
nozzle actuation system is established; The load spectrum of the actuating system under different working 
conditions and deflection angles is obtained by CFD numerical simulation results; Using AMESim, the 
hydraulic actuation system model of two-dimensional vector nozzle is established. Finally, a two- 
dimensional vector nozzle/actuator/turbofan engine comprehensive model is established based on the 
component-level model of turbofan engine. The comprehensive model can accurately calculate the changes 
of main performance parameters such as engine rotor speed, thrust and pre-turbine temperature when the 
nozzle vector deflects. The simulation results show that the load on the nozzle actuator is positively 
correlated with the deflection angle, and the load change will lead to the fluctuation of the output 
displacement of the hydraulic actuation system. During vector deflection, due to the change of engine 
circulation capacity and thrust coefficient, the performance parameters such as engine low-pressure rotor 
speed and surge margin fluctuate greatly. 


INDEX TERMS 2DCD, Actuator load, Aircraft engine mode, Comprehensive model 


I. INTRODUCTION 

With the development and application of thrust vector 
technology, the flight performance, take-off and landing 
performance, maneuverability and post-stall maneuverability 
of aircraft have been greatly improved, so advanced aviation 
countries in the world give priority to it as an important 
technology [1]. Vector nozzle technology is the foundation 
and core of thrust vector technology. Among various vector 
nozzle schemes, 2D vector nozzle (2DCD) has the 
characteristics of low infrared radiation and easy integration 


two-dimensional nozzle design and ground tests, and on this 
basis, carried out ground and flight verification tests of 
flying-weight nozzles [5]. Since 1989, the two-dimensional 
vector nozzle has been tested on the F-15. In the 1990s, 
thrust vector aircraft made great progress in engineering, and 
a new type of vector thrust aircraft was studied [6]. After 
long-term theoretical and experimental research of thrust 
vector technology and the accumulation of several 
demonstration test projects, at the beginning of the 21st 
century, the thrust vector technology was successfully 


design with rear fuselage [2], so 2D vector nozzle is an 
essential technology for future advanced fighters [3]. 

Since the 1970s, NASA, military, aircraft and engine 
companies in the United States have carried out a series of 
two-dimensional vector nozzle research plans and wind 
tunnel test plans for thrust vector aircraft [4]. In the 1980s, 
the United States carried out a full-scale, non-flying-weight 
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applied to the main fighters and stealth aircraft of the Air 
Force. In the academic field, scholars at home and abroad 
have done a lot of work in numerical simulation and 
experiments on the internal flow characteristics of two- 
dimensional vector nozzle [7-15]. However, the research on 
actuator load and vector nozzle/engine integrated model is 
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mostly focused on axisymmetric vector nozzle. Zhang [16] 
put forward a coupling modeling method of axisymmetric 
vector nozzle spatial kinematics model and vector deflection 
CFD model, which solved the problem that the load of A9 
actuator is unknown when axisymmetric vector nozzle 
actually works. Huang [17] systematically modeled the 
vector nozzle control device based on AMESim, and studied 
the steady-state performance, dynamic performance, 
robustness and performance when the fault returns. As a 
necessary component of the engine, the working state of the 
engine fluctuates violently when the vector is deflected, so 
Cai [18] have established a dynamic model of turbofan 
engine with axisymmetric vector nozzle, and the linear active 
disturbance rejection control is adopted to reduce the 
variation amplitude of engine state parameters under the 
traditional PID closed-loop control. Qu studied the influence 
of vector deflection on the working state of the engine under 
open-loop control; Li [20] established a micro turbojet vector 
propulsion system model based on experimental data and 
mechanism model. Li's model based on experimental data is 
difficult to ensure the accuracy requirements, Cai and Qu 
only established the mathematical model of the engine with 
axisymmetric vector nozzle. However, the structure and 
kinematics model of two-dimensional vector nozzle are quite 
different from that of axisymmetric vector nozzle. It is 
impossible to obtain the influence of the vector deflection of 
the two-dimensional vector nozzle on the working state of 
the engine. Therefore, it is of great significance to establish a 
comprehensive model of two-dimensional vector 
nozzle/actuator load/turbofan engine based on mechanism for 
the stability and accurate control of engine and nozzle in the 
next step. 

In this paper, the aerodynamic characteristics of the two- 
dimensional vector nozzle are combined with the mechanical 


Ag actuator 


motion mechanism of the nozzle, and the internal flow 
characteristics of the two-dimensional vector nozzle under 
different working conditions and different deflection angles 
are calculated by CFD. Based on the spatial kinematics 
model of the two-dimensional vector nozzle, the mechanical 
model of the vector nozzle actuation system is established by 
combining force analysis and torque analysis, and the 
hydraulic actuation system model of the nozzle is established 
by AMESim, and a comprehensive model is established by 
combining with the component-level model of the turbofan 
engine. 

Il. Mechanical model 

The structural diagram of the two-dimensional vector nozzle 
is shown in Figure 1, and its convergent section and 
divergent section can be continuously adjusted. The upper 
expansion adjusting plate and the lower expansion adjusting 
plate are respectively driven by an A9 actuator, the front end 
of the A9 actuator is hinged with the casing, and the rear end 
is hinged with the square ring connecting rod mechanism. 
The three-dimensional structural diagram of the square ring 
linkage mechanism is shown in Figure 2, and the square ring 
can only move horizontally in the groove on the side wall of 
the casing, thus ensuring that the A9 actuator always 
maintains a horizontal position. The rear end of the square 
ring linkage mechanism is hinged with the expansion 
adjusting plate. When the output displacement of the A9 
actuator changes, the square ring moves horizontally, thus 
driving the pull rod and the expansion adjusting plate to 
move. When the upper and lower A9 actuators 
synchronously output displacement, adjust the A9 cross- 
sectional area; When the output displacements of the upper 
and lower A9 actuators are different, the section of A9 tilts to 
realize pitching thrust. 


Square ring » 


OV Q 


Round to square 
transition section 


7 section 


8 section 9 section 


FIGURE 1. Schematic diagram of spatial motion and force of two-dimensional vector nozzle 
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A9 actuator 


FIGURE 2. Structural sketch of expansion adjusting plate 


The A8 actuating system of two-dimensional nozzle is 
connected with A8 actuating cylinder and convergence 
adjusting plate by rocker arm, and its structure diagram is 
shown in Figure 3. The front end of the rocker arm is hinged 
with the A8 actuating cylinder, and the rear end is hinged 


FIGURE 3. Structure diagram of convergence adjusting piece 


The force acting on the actuating system of two- 
dimensional vector nozzle is shown in Figure 1, and the side 
plate of the two-dimensional nozzle is fixed and cannot be 
adjusted. When analyzing the force acting on the actuating 
system, the friction between the upper and lower adjustable 
plates and the side plates will not be considered for the time 
being. The force acting on the expansion plate includes the 
aerodynamic resultant force F of the expansion plate, the 
force Fag exerted on it by the pull rod AB and the force Fr 
exerted on it by the convergence plate. 

When the braking system is in a steady state, it can be 
obtained from the moment balance between the convergent 
adjusting plate and the divergent adjusting plate at point D: 


F,,xd=FxX, (1) 
d is the distance from the pull rod AB to the point D, and Xc 
is the position of the aerodynamic resultant force in the 


expansion section. 
d can be calculated by equation (2): 
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Upper expansion 
regulating sheet 


with the convergence adjusting piece through the pull rod. 
When the A8 actuator outputs displacement, the rocker arm 
rotates around its rotating rod, thus driving the pull rod and 
the convergence adjusting plate to rotate, thus realizing the 
change of A8 throat area. 


ee 


ea ee ore sae regulator 


+ 
tan(v + ô) 


d =sin(v + 0): i be (2) 


v is the expansion angle of the expansion adjusting plate, ô is 

the included angle between the pull rod AB and the 

horizontal axis, and L; is the distance from point B to point C. 
6 calculation equation is: 


(3) 


5- arcsin Ž: —R, — Lpe ‘sin v — L, se 


L 


‘AB 


Ry is the distance from point A to the nozzle axis, Rp is half 
the throat height, Lps is the distance from point D to point B, 
and Laz is the length of pull rod AB. 

The force component of the pull rod AB along the 
horizontal direction is balanced by the A9 actuator, and then 
the output force of the A9 actuator can be solved as follows: 


Fy = Fy, cos (4) 


The acting force Fr of the expansion adjusting plate on the 
convergence adjusting plate is obtained by Equation (5): 
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F, =-(F + Fig) (5) 


To simplify the calculation, decompose it along the X and 
Y axes, and we can get: 


F, =—(-F -sinv + F,, -cosd) © 
F, = -(F -cosv — F,,-sind) 

The force exerted by the pull rod GH on the convergence 
adjusting plate can be obtained from the moment balance of 
the connection point E between the convergence adjusting 
plate and the casing: 


FIGURE 4. Enlarged view of convergent section structure 


When the actuating system is in a stable state, the moment 
balance between A8 actuating cylinder and connecting rod 
GH about fixed point N can be obtained as follows: 


Ferd Eyxd (8) 


Fag is the output force of actuator A8, dz is the distance from 
actuator MK to fixed point N, and d; is the distance from rod 
GH to fixed point N. 

də can be expressed by Equation (9): 


d, = Lpy ‘sin 0 (9) 


Where Lry is the length of the left KN of the rocker arm, and 
0 is the included angle between the actuator MK and the left 
KN of the rocker arm. 
0 can be represented by equation (10): 
Lx + Lx z Lan 
10 
=m) (10) 


‘NK 


0 =arccos( 


` MK 


Where Lyn is the distance between point m and fixed point n 
of rocker arm. 

The distance of the dotted line NG can be expressed by 
Equation (11): 


Lyc = Lro + Lin — 2+ Lpg + Ley cose + G) d1) 


Where Leg is the distance between point E and hinge point G, 
Ley is the distance between point E and fixed point N, and ¢ 
is the included angle between connecting line EN and the 
central axis. 


Ag actuator 


Fp + Lgy sine + Fy Lep cose + P+ Xp 
Fon = d 


(7) 


1 


dı is the distance from the point E to the GH point, Lep is 
the length of the convergence adjusting plate, ¢ is the 
convergence angle of the convergence adjusting plate, and Pz 
and Xz are the aerodynamic resultant force and center of the 
convergence adjusting plate. 

The output force of A8 actuator can be further solved by 
the force acting on the pull rod GH. 


PA Rotating lever 
: © 


Any 


The included angle y between the right end NH of the 
rocker arm and the rod GH can be expressed by Equation 
(12): 

2 2 _ 72 
y= arccos( e Eeu Pas + Lon -Lwa ) (12) 
‘Lyn Lon 
Where Lac is the length of the rod NG and Lug is the length 
of the right half of the rocker arm. 

Then the distance d; from the rod GH to the fixed point n 

can be expressed by equation (13). 


d, =Lyg siny (13) 


The expression of the output force F'4s of A8 actuator can 
be obtained by bringing equations (9) to (13) into equation 
(8). 

The included angle J between the rod HG and the dotted 
line NG can be expressed by Equation (14): 


‘NG 


£ =arccos( 
“SHG 


The included angle a between rod EG and dotted line NG can 
be expressed by equation (15): 


a = arcsin -sin(£ + €)) (15) 


‘GN 


Then the distance d1 from the HG point to the E point of the 
pole can be expressed by equation (16): 


d, = Leg ` sin(a + p) (16) 
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By introducing the equation (7), the force exerted by the pull 
rod GH on the convergence regulator can be obtained. 


Ill. Actuator load model 
The two-dimensional vector nozzle consists of a circular 
torque transition section and a main nozzle section 
(adjustable convergence section and adjustable expansion 
section), and its simplified model is shown in Figure 5. The 
circular torque transition section is not adjustable, which 
realizes the transition from the circular outlet of afterburner 
to the rectangular inlet of the main nozzle section. The main 
nozzle section consists of two parallel side plates, left and 
right, and two adjustable adjusting plates, upper and lower. 
Both convergent section and divergent section can be 
adjusted to realize the change of throat area and outlet area of 
two-dimensional vector nozzle, so as to meet the 
requirements of aircraft under different working conditions 
and provide vector thrust and reverse thrust. H7, Hs and Ho 
are the heights of the inlet, throat and outlet of the main 
nozzle respectively; W is the width between the two side 
walls; a is the convergent half angle; J is the expansion half 
angle; R,, is the radius of throat arc; L, represents the length 
of the circular torque transition section [21]. 
Circular totini Adjustable 
transition Adjustable extension 


section cohvergence segment 
egment 


Narrow-sided 
symmetrical 

plane section 
DY 


W/2 

Wide-sided 
symmetrical 
lane section 


FIGURE 5. 2DCD 2D schematic diagram 

The grid model of two-dimensional vector nozzle is shown 
in Figure 6. In order to reduce the difficulty of grid division 
and ensure the accuracy of calculation, unstructured grid is 
used in the transition section of circular torque, structured 
grid is used in the convergence and expansion section, and 
encryption is carried out near the nozzle wall. 


FIGURE 6. Grid division 

CFD numerical simulation of vector deflection of two- 
dimensional vector nozzle is based on Fluent software, and 
S-A single equation model is selected as turbulence model 
[22]. The nozzle inlet is given the pressure inlet boundary 
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condition, and the nozzle outlet is given the pressure outlet 
inlet. The specific values are determined by the engine model. 
The nozzle thrust calculation equation is as follows: 


F, = Amu+(p- p, )AA, 
F, =$ Amv+(p- p,)^, (17) 
F, = Amw+(p- p,) AA. 


Where F», Fy and F, respectively represent axial thrust, yaw 
thrust and pitch thrust; Am is the mass flow of each micro- 
element area of the outlet section A9. 

The effective vector angle fo is: 


Ay =arctan( F/F.) (18) 


The total thrust F is: 


F-J (19) 


The thrust coefficient Cm is the ratio of the actual total 
thrust to the ideal total thrust at the nozzle outlet section, 
namely: 


C,=F/F, (20) 
The aerodynamic resultant force of the expansion tab is: 
F= j” Pœ): Sœ)dx (21) 


Where P(x) is the static pressure distribution along the axial 
center line of the expansion tab, and the fitting curve of the 
static pressure distribution at the center line of the upper and 
lower expansion tabs of the two-dimensional nozzle in each 
deflection state obtained by Fluent calculation is solved, 
where S(x) is the distribution along the axial area of the 
expansion tab and the integral interval is the length of the 
expansion section. 

The position of the action point of aerodynamic resultant 
force in the expansion section is: 


1 Lpp 
m=z Í, P(x) + S(x) + xdx (22) 


The aerodynamic force of that convergence regulator are 
as follows: 


P= j” P(x)-S(xddx (23) 


In the equation, the upper integral limit Lep is the length 
of the convergence regulator, and P.(x) is the aerodynamic 
load distribution along the convergence regulator. Similarly, 
it can be calculated and solved by vector deflection Fluent. 


IV. Comprehensive model 

Hydraulic servo control system is a feedback control system 
with hydraulic power components as driving devices [23,24]. 
In this system, the output (displacement, speed, force, etc.) 
can automatically, quickly and accurately reproduce the 
changing law of the input, and it is also a power amplifier 
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[25,26]. In this section, based on the load characteristics, the 
simplified model of A8 actuating system including 
electromagnetic directional valve, hydraulic actuator, 
displacement sensor, external load force model and PID 
controller is established by using AMESim software. 


FIGURE 7. A8 actuation system model 

Set the input parameters of A8 actuating system as 
displacement command and external load force, and the 
output parameters as displacement feedback from 
displacement sensor. The established AMESim model is 


| Supersonic Inlet 


Compressor 


22 


shown in Figure 7. Set the input parameters of A8 actuating 
system as displacement command and external load force, 
and the output parameters as displacement feedback from 
displacement sensor. The established AMESim model is 
shown in Figure 7. 

The engine component-level model in this paper is a twin- 
shaft mixed-row turbofan engine with small bypass ratio. 
Figure 8 shows the overall structure of the engine, and its 
main components include supersonic inlet, fan, compressor, 
external bypass, combustion chamber, high-pressure turbine, 
low-pressure turbine, mixing chamber, afterburner and tail 
nozzle. 2~22 sections are fans, 13~16 sections are external 
ducts, 25~3 sections are compressors, 3~41 sections are 
combustion chambers, 41~45 sections are high-pressure 
turbines, 45~5 sections are low-pressure turbines, 6~7 
sections are mixing chambers and afterburner, 7~9 sections 
are tail nozzles, and 8 sections are nozzle throat areas. 


Bypass H-turbine 


35 rá Nozzle 
13) 3/41 16 6 65 
an i 


Inlet Ramp 
FIGURE 8. Schematic diagram of turbofan engine structure 


The component-level model of aero-engine dynamic 
process includes component model and solution of common 
working equation. The component model simulates the 
aerodynamic and thermodynamic process in the component, 
which is not described in detail in this paper. The matching 
relationship between components in the steady-state process 
is determined by the relationship of continuous flow, power 
balance and pressure balance. In this paper, the relative 
physical speed Nz of the low-pressure rotor, the relative 
physical speed Ny of the high-pressure rotor, the fan pressure 
ratio coefficient Z; the high-pressure compressor pressure 
ratio coefficient Zey, the high-pressure turbine pressure ratio 
coefficient Z; and the low-pressure turbine pressure ratio 
coefficient Z are selected as unknown parameters to solve 
the component-level model, and six unknown parameters 
correspond to six independent equations, as shown in 


Equation (24). 
é(U)=f(N)=0 (i=1~6) (24) 


Where N represents a vector consisting of six initial guesses, 
which can be expressed as 


N =[n,n,n,n, nn] 
7 (25) 
|N; N 2 Zar Za Ziy | 
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The six equilibrium equations are as follows: 

(1) The residual equation of the balance between the fan 
consumption power Wp and the power Wi, provided by the 
low-pressure turbine is obtained from the power balance of 
the low-pressure rotor: 

Wit. — W, 


ga- t S 26 
' Wan = 


Where 7 is the efficiency of the low-pressure rotor. 

(2) Based on the power balance of the high-pressure rotor, 
the residual equation of the balance between the power 
consumption W.y of the high-pressure compressor and the 
power output Wew of the high-pressure shaft and the power 
W, provided by the high-pressure turbine is obtained: 


az Waa — Wor + Woa) 
í Waa 


(27) 


Where 7 is the efficiency of the low-pressure rotor. 

(3) Based on the flow balance of high-pressure turbine, the 
actual inlet flow m4, of high-pressure turbine is obtained and 
the residual equation of turbine flow m4cx is obtained by 
three-dimensional interpolation with the characteristics of 
high-pressure turbine: 
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& = My. = Ng) ox (28) 


Maie 


(4) Based on the low-pressure turbine flow balance, the 
actual low-pressure turbine inlet flow mas, and the low- 
pressure turbine characteristics are interpolated to obtain the 
residual equation of turbine flow m4s5,cx: 

_ Mas, — Mas ex 


E, = — Mma (29) 


M45 


(5) Based on the static pressure balance at the outlet of the 
external culvert, the residual equation between the internal 
static pressure Ps of the engine and the static pressure P36 of 
the external culvert is obtained: 


_ ee 
P 


6 


(30) 


és 


(6) The residual equation of the nozzle throat static pressure 
Ps and the static pressures Ps. calculated according to the 
internal flow characteristics of the nozzle is obtained from 
the nozzle throat static pressure balance: 

RK. 


a (31) 


8 


In the dynamic process of the engine, the rotor power of 
the high and low pressure turbine is no longer balanced with 
the rotor power of the compression part, so the deviation 
value will change the rotor speed, and the balance equation is 
replaced by the dynamic equation of the high and low 
pressure rotor. 


dN zy 
TF = (Wani ~ voll m (=| | 


dN zY 
P = (Wahu — Wor -Wa Nery (=) | 


(32) 


Where J; and Jy are the moment of inertia of the low 
pressure shaft and the moment of inertia of the high pressure 
shaft respectively, and their values are obtained with 
reference to Gasturb software. In the dynamic process, the 
current rotational speed (N2) is obtained by superposition of 
the previous rotational speed (V7) by Euler method, namely: 


ie =N,,+dN,At (33) 


Np = Np, + dN, At 


In order to obtain the aerodynamic thermodynamic 
parameters and performance parameters of each section of 
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the engine, the Newton-Raphson method with fast 
convergence speed and small amount of calculation is used to 
iteratively calculate the balance equation. 

The effect of vector deflection on engine performance is 
mainly reflected in the change of effective throat area and 
nozzle thrust coefficient. The vector deviation transfer causes 
the nozzle throat to tilt, resulting in the reduction of the 
effective throat area, thus affecting the engine working state, 
while the change of thrust coefficient causes the change of 
the engine output thrust. Therefore, the effective throat area 
coefficient Case is defined as: 

C, = (84) 


Ase p 
Maxi is 


Where m, is the nozzle flow calculated by CFD in the 
undeflected state and m, is the nozzle flow calculated by 
CFD in the deflected state. 

Therefore, the calculation equation of Ps: in the static 
pressure balance equation of the nozzle throat is as follows: 


Ms lT, > 
P, g8 84 ( Ag ) (35) 
Define relative thrust coefficient Cmr: 
Car = Gin! Garis (36) 


Where Cm is that thrust coefficient of the deflect nozzle 
under different work conditions, and Cypaxis is the reference 
thrust coefficient of the undeflected nozzle. 

When the vector deflects, the equation for calculating the 
resultant thrust of the engine vector is: 


F = (ting “Vg — Myo “Vo + (Fs P): Ay) Cpa (37) 


When the vector deflects, the effective throat area 
coefficient Cy4se changes, which changes the common 
working equation, while the relative thrust coefficient Cmr 
changes, which leads to the decrease of the thrust of the tail 
nozzle, and both of them work together to affect the change 
of the engine performance. The control plan of controlling 
the high-pressure rotor speed with the main fuel and the low- 
pressure rotor speed with the nozzle throat area A8 is adopted, 
and the vector nozzle actuator model is combined to establish 
a two-dimensional vector nozzle/engine comprehensive 
model as shown in the structural block diagram of Figure 9. 
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FIGURE 9. Structure block diagram of vector deflection comprehensive model 


V. Simulation and Analysis of Vector Deflection 


A. Aerodynamic performance simulation of two- 
dimensional vector nozzle 
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The simulation results of aerodynamic performance of two- 
dimensional vector nozzle under different working 
conditions and different vector deflection angles are shown in 
Figure. 
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FIGURE 10. Structure block diagram of vector deflection comprehensive model 


As can be seen from Figure 10, with the increase of the 


vector deflection angle, the engine flow coefficient decreases. 


When the deflection angle is less than 10°, the vector 


deflection has little influence on the nozzle flow performance. 


When the deflection angle is greater than 10°, the nozzle 
flow performance decreases to a certain extent, and the 
changing trend is basically the same in different flight States. 
When the vector deflection angle is 20°, the nozzle flow 
capacity decreases by about 6.5%, and the engine state is 
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greater at the same deflection angle. the effective vector 
angle of the engine is slightly smaller than the vector 
deflection angle, which is especially obvious in the high- 
altitude flight state. Similarly, the greater the engine state, the 
closer the effective vector angle is to the vector deflection 
angle, and the better the linearity between them. 

with the increase of the vector deflection angle, the total 
pressure recovery coefficient of the engine decreases. When 
the deflection angle is less than 5°, the vector deflection has 
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little influence on the total pressure loss of the nozzle, and 
when the deflection angle is greater than 5°, the total pressure 
loss decreases to a certain extent, and the change trend is 
basically the same in different flight states. Vector deflection 
has little effect on thrust coefficient. Under different flight 
conditions and different engine conditions, with the increase 
of vector deflection angle, thrust coefficient shows a slight 
downward trend. When the vector deflection angle is 20°, the 
relative thrust coefficient of the nozzle decreases by about 
0.5%. 


B. Simulation of load characteristics of actuator 
According to the principle of hydrodynamics, the pressure on 
the wall in the flowing pipeline is equal to the pressure at the 
point where the velocity is 0. Therefore, it is necessary to 
create a central surface in the process of extracting the wall 
pressure, extract the point on the central surface with the 
velocity of 0 for wall pressure function fitting, and import the 
upper and lower wall data in supersonic cruise state into the 
fitting software. The fitting result of the upper wall pressure 
is shown in Figure 11 and 12. 

In order to calculate the load characteristics of the 
actuating system under different working conditions and 
deflection angles, this paper selects the ground point 
PLA=70°. 90°. 110°, the high altitude point H=10km, 
Ma=0.8, PLA=70° subsonic cruise state, and the supersonic 
cruise state H=10km, Ma=1.2, PLA=110°. The vector 
deflection angles are 0°. 5°. 10°. 15° and 20° 
respectively. Figures 12 show the nozzle convergence and 
the loading conditions of the expansion section regulator 
under different working conditions, Figure 13 shows the 
loading conditions of A8 and A9 actuators. 
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FIGURE 14. Loading conditions of actuators A8 and A9 in different states and vector deflection angles 


In the process of vector deflection, the resultant force on 
the upper and lower clamping plates in the convergence 
section is basically unchanged. This is because the air flow 
accelerates the expansion in the binary vector nozzle and 
reaches the sound velocity at the throat. The impact of the 
vector deflection in the expansion section will not spread 
forward. The resultant force on the upper adjusting plate 
increases continuously, while the resultant force on the lower 
adjusting plate decreases slightly. 

With the increase of the vector deflection angle, the load 
force of the upper A8 actuator increases continuously, and 
the load force of the lower A8 actuator first increases and 
then decreases. The load force of the upper A9 actuator 
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decreases and the load force of the lower A9 actuator 
increases. 


C. Hydraulic actuator simulation 

According to the calculation results of load force of A8 
actuator in the state of vector deflection in the previous 
section, the vector deflection simulation of A8 actuator 
system is studied. First, input the command of centering after 
two 5cm steps as shown in Figure 15 and the command of 
constant A8 actuator load 50KN as shown in Figure 14, and 
the displacement response result of the actuator is shown in 
Figure 14. 
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FIGURE 15. Simulation of vector deflection of actuating system 


From fig. 15, it can be concluded that the small closed-loop 
control loop of the hydraulic actuation system of the binary 
vector nozzle actuator established by AMESim software 
based on the load analysis results has fast response and no 
overshoot. According to the response results, with the 
increase of the vector deflection angle, the load force of the 
A8 actuator changes to some extent, which will lead to 
certain displacement fluctuation of the A8 actuator under the 
control of the constant displacement command. The larger 
the vector deflection angle, the more drastic the change of the 
load force of the A8 actuator, resulting in the greater the 
change amplitude of the output displacement of the actuator, 
and the displacement change of the A8 actuator at the high 
altitude point is smaller than that at the ground point. 


D. Integrated model simulation 

Open-loop simulation of the established engine component- 
level model is carried out, and the instruction of ramp fuel 
supply law as shown in Figure 17 is input. The relative 
physical speed and thrust response of the low-pressure rotor 
of the engine basically meet the expected dynamic and 
steady-state performance of the engine. 

In the high-altitude cruise state, the vector deflection 
signal shown in fig. 17(a) that quickly deflects 20 degrees 
and then returns to the center is input to the comprehensive 
model, and the load change instruction corresponding to the 


vector deflection shown in fig. 17(a) is input at the same time. 


The simulation results are shown in figs. 17 (b)~(h). 
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FIGURE 16. Open loop simulation 

From Figure 17(b)~(f), it can be seen that the response of 
engine performance parameters is consistent with the 
theoretical calculation and aerodynamic simulation results. 
When the vector deflection increases, the effective throat 
area decreases and the engine flow decreases, which leads to 
the rapid decline of engine thrust and speed. Under the action 
of closed-loop controller, the engine fuel supply rate 
increases and the nozzle expands to increase the inlet flow. 

As can be seen from Figure 17(g)~(h), vector deflection 
has little influence on compressor rotor speed and turbine 
front temperature, but the response of engine speed and other 
parameters fluctuates greatly, with the maximum variation of 
fan speed of about 1.78% and high pressure surge margin of 
about 2.88%. 
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VI. CONCLUSION 

In this paper, a two-dimensional vector nozzle/actuator 
load/turbofan engine comprehensive model is proposed and 
established, and the main conclusions are as follows: 

(1) During vector deflection, the resultant force on the upper 
and lower tabs of the convergent section of the two- 
dimensional vector nozzle remains basically unchanged, 
while the resultant force on the upper tab of the divergent 
section increases continuously, while the resultant force on 
the lower tab decreases slightly. 

(2) With the increase of vector deflection angle, the load 
force on A8 and A9 actuators fluctuates greatly, among 
which the load force on the upper A8 actuator changes the 
most. When the vector deflection is 20, the load force on the 
upper A8 actuator increases by about 2 times, and the output 
displacement fluctuates by 1.225%, which affects the stable 
control of the vector deflection process A8. 

(3) Vector deflection will change the effective throat area and 
thrust coefficient of the engine nozzle, which will lead to 
severe fluctuation of engine performance parameters. The 
maximum variation of fan speed is about 1.78%, and the 
maximum variation of high pressure surge margin is about 
2.88%, but it has little effect on compressor rotor speed and 
turbine front temperature. 
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